Manifestation of anisotropy persistence in the hierarchies of MHD scaling exponents 
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The first example of a turbulent system where the failure of the hypothesis of small-scale isotropy 
restoration is detectable both in the 'flattening' of the inertial-range scaling exponent hierarchy, and 
in the behavior of odd-order dimensionless ratios, e.g., skewness and hyperskewness, is presented. 
Specifically, within the kinematic approximation in magnetohydrodynamical turbulence, we show 
that for compressible flows, the isotropic contribution to the scaling of magnetic correlation functions 
and the first anisotropic ones may become practically indistinguishable. Moreover, skewness factor 
now diverges as the Peclet number goes to infinity, a further indication of small-scale anisotropy. 
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A wide interest has recently been devoted to the pos- 
sible occurrence of small-scale isotropy restoration for 
scalar (see e.g., Ref. [Q and references therein), Navier- 
Stokes and magnetohydrodynamical (MHD) turbu- 
lence pppL The scenario can be summarized as follows. 
In the presence of anisotropic large-scale injection mech- 
anisms, the inertial-range statistics is characterized by 
an infinite hierarchy of scaling exponents; however the 
leading contribution to scaling comes from the isotropic 
component. From this point of view, one might argue 
that large-scale anisotropy does not affect inertial-range 
scaling properties. Actually, focussing on a larger set of 
observables, small-scale anisotropies become manifest. It 
turns out that the behavior of odd-order dimensionless 
ratios (e.g., skewness and hyperskewness) is completely 
different from the case of small-scale isotropy restora- 
tion. Such indicators go to zero down to the inertial 
range much slower than predicted by dimensional con- 
siderations or, more dramatically, they diverge at the 
smallest scales || (see also [|]||). 

The main aim of this Letter is to present a model of 
MHD turbulence where, by varying the degree of com- 
pressibility of the velocity field, anisotropic persistence 
is now detectable both from the 'flattening' of the hier- 
archy of inertial-range scaling exponents (the isotropic 
component and the first anisotropic ones may become 
practically indistinguishable) and the divergence of skew- 
ness factor with the Peclet number. We give the basic 
ideas and results; longer and more exhaustive technical 
discussions will be presented in a forthcoming paper. 

In the presence of a mean component B° (actually 
varying on a very large scale L) and for the compressible 
velocity field v, the kinematic MHD equations describing 
the evolution of the fluctuating (divergence-free) part B 
of the magnetic field are jnj: 



d t B a +v-8B a = -(B a + B° a )d ■ v + B ■ dv a 



where kq is the magnetic diffusivity. The field B° plays 
the same role as an external forcing driving the system 
and being also a source of anisotropy for the magnetic 
field statistics. 

Our choice for the velocity statistics generalizes that of 
the well-known kinematic Kazantsev-Kraichnan model 
for the compressible case: v is a Gaussian process of zero 
average, homogeneous, isotropic and white in time. It 
is self-similar and defined by the two-point correlation 
function: 

{v a {t 1 x)vp{t',x'))=5(t-t')[d%-S a p{x-x% (2) 

where d^g = const 6 a g and S a g(x — x') is fixed by 
isotropy and scaling: 
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The degree of compressibility is thus controlled by the 
ratio p = C 2 /S 2 , with S 2 oc ((diVkdiVk)) and C 2 cx 
((d ■ v) 2 ). It satisfies the inequality < p < 1; p = 
and 1 corresponding to the purely solenoidal and poten- 
tial velocity fields, respectively. 

In the present Letter our attention will be focused on 
the inertial-range behavior of magnetic correlation func- 
tions, where power laws arc expected in their decompo- 
sitions on a set of orthonormal functions P 



(Brf(t,x)Bf,(t,x')) =E£o-Pi(«*< 
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Bu being some component of B, e.g., its projection along 
the direction r = r/r or B = B°/B°, and <p is the angle 
between r and B°. Notice that, due to the anisotropic 



injection mechanism, the inertial-range statistics is now 
characterized by an infinite hierarchy of exponents (j de- 
notes the j-th anisotropic sector), rather than just one 
exponent as in the isotropic case. 

If the Kolmogorov 1941 isotropization hypothesis 
holds, the contribution from the anisotropic sectors 
(i.e. j ^ ) to the scaling of correlation functions should 
be negligible with respect to the isotropic component. 
Such picture indeed holds for even correlation functions 
and solenoidal velocity ||]6| , but it breaks down for com- 
pressibility strong enough. 

In order to prove this fact, consider first the pair cor- 
relation function, C a p(t,r) = (B a (t, x)Bg(t, x')). Here, 
exploiting the zero- mode technique 11 - 14 1, the complete 



set of scaling exponents q ' can be found nonperturba- 
tively. We give the basic ideas of the strategy. 

The first key point is that exact, closed equation for 
C a p can be found due to the time decorrelation of the 
velocity field: 



dt C a p — Sij didj C a p — (dj Sip) di C a j 

{djS ai ) {diCjp) + (d^Sap) (Gy + 



2K d 2 c a p - BZBpidjSd - irir;<),<).s,, + 

B a Bpdidj &ij ~ C a j {didj Spf) + Cpj didj S a p 
C a p didj + 2(diC a p) dj Sy 



(6) 



(see H for the derivation in the incompressible case). 
C a p can be projected on a vector basis including the 
large scale magnetic field B° [f|](|, or, alternatively, on 
the purely inertial range basis that span the irreducible 
representations of SO(g?) [fl5|| . Deep in the inertial range, 
both descriptions yield a system of linear algebraic equa- 
tions for the scaling law coefficients. The leading solu- 
tions are associated with the homogeneous solutions of 
such system, i.e. with zero modes. Their scaling expo- 
nents follow from the imposition that the determinant of 
the coefficients be zero. Schematically, the solutions can 
be expressed in the form: 



Cj' 1 = a + \J j3 + -fVS (j even), 



(7) 



where a, 0, 7 and S are cumbersome functions of £, d and 
j and will not be reported here for the sake of brevity. 
Contributions to the scaling associated with odd j's van- 
ish due to the symmetry <j) 1— * — <j). 

In particular, the expression derived in Ref. |p^| in the 

isotropic case has been recovered here for j ' = 0. 

For j = and j — 2, the limit £ — > in Eq. (Q) yields: 



c ia = (-1 + 2 p - d p)z+o(e), 
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In Fig. |l| we present the behavior of obtained from 
the expression (]?]) for p = (thin lines) and p = 1 
(heavy lines). As one can see, for p = we have 
Co' 1 < but C2' 1 > and C4' 1 > so that, in partic- 
ular, (r/Lpa ^> (r/Lps in the inertial range of scales. 




FIG. 1. Behavior of q' vs £ for d = 3 and j = 0, 2 and 4 
(from below to above). Heavy lines: p = 1; thin lines: p = 0. 

The situation changes for p > p c (p c ~ 0.1274 for 
j = 2 and d = 3), when £ 2 ' becomes negative for all 
< £ < 1. This means that, for compressibility strong 
enough, (r^Lp° and (r/Lp* become very close (still, 

< Ck' 1 f° r 3 < As we can see f rom Eq s - (|ic|), 
(O) the effect becomes dramatic for large dimensions and 
p ^ when C^' 1 ~ C2' ■ The contribution to the scaling 
of the correlation functions in Eq. (|5|) coming from the 
sector j = 2 thus becomes less and less subleading as p 
and/or d increase. 

Further strong evidence of the crucial role of compress- 
ibility in the failure of the small-scale isotropy restoration 
can be obtained by looking at the higher-order correlation 
functions [i.e. n and/or q > 1 in Eq. (^|)] and, in particu- 
lar, at dimensionlcss ratios of odd-order moments. Were 
isotropy restored at small scales, such ratios would go to 
zero for large Peclet number. The latter is defined as Pe 



(L/rj) 1 ^, where L is the integral scale and r] 



cx K, 



while for large d we have: 



the dissipation scale. 

Non-zero values of such indicators are thus the signa- 
ture of anisotropy persistence. As we are going to show, 
in contrast to the incompressible case, for p large enough 
the skewness factor now diverges as Pe — > 00. To be more 
specific, the leading contribution in the expression (^) can 
be written as: 

(If (t, x)B*(t, *')> cx (r/r,) a o 9 (r/L)^' q cx r-C 9 , (12) 

where we have defined £q ' q = «g' 9 + /?q' 9 . The expres- 
sions for £q ' 9 have been obtained up to first order in £ by 
means of the field theoretic renormalization group and 
operator product expansion. A detailed presentation of 
these techniques for the case p = can be found in 



(see also Refs. |||]|l7]] for the scalar case and Q for 
general review); below we confine ourselves to only the 
necessary information. 

The key role is played by the scaling dimen- 



sions A[n,j] of the j-th rank composite operators 
B ai (x) ■ ■ ■ B aj (x) [B a (x)B a (x)] ' (n = 21 + j is the total 
number of B's). The first order in £ (one-loop approxi- 
mation) yields 



A[n, j] = " {n(n - 1) [2 - p(d 3 ~d 2 -2d + 4)] - (n - j)(n + j + d - 2)(d + 1 - 2p)} + 0(f)- (13) 



The exponents in Eq. (c 
( |l3| ) as follows (see Ref. | 



are related to the dimensions 



A[n + q,j]-A[n,j n ]-A[q,j g ], 



where j n = for n even and 1 for n odd. We thus 
conclude that the inequality dA[n, j]/dj > 0, which 
follows from Eq. (13) for all p and d > 2, general- 
izes the hierarchy discussed above to the higher-order 
functions. It becomes flatter and flatter as p grows 
(d 2 A[n,j]/djdp < 0), while for d — > oo (and p ^ 0) 
the effect becomes even stronger: in the leading 0(d), 
expressions ( p"3"| ) and ( |l4| ) are now independent of j. 

The leading term of Eq. (^|) in the inertial range is 
given by the contribution with the minimal j: 



-i 



nq(l + pd 2 - 2p) + (d+l - 2p) 



(d + 2\ 



(d + 2) 
■(l + pd 2 -2p), 



(14) 



where the first holds if both n and q are odd, 
and the second otherwise. For n = q = 1, 
expression (|J) is recovered. Knowing the expo- 



nents a, 



n . q 



and /3, 
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dimcnsionless ratios of the form 



R 



2n+l 



(r) ee (B 2n (x)B ll (x')}/(B ll (x)B ll (x')}( 2n+i y 2 



can 



be constructed and, as in Ref. M, e valuated at the dis- 
sipative scale [i.e. r = rj in Eq. (|12[)]. When doing this, 
the explicit dependence on Pe appears and the final 0(£) 
expressions read: 



i? 2 „+i(Pe)cxPe' 
with the exponents 



C2n + 1 



1,2,-' 



(15) 



^_ 2n 2 (l-2p + pd 2 ) (l-2p+pd) 
a 2n+1 (p) = ^ . (16) 

It is easy to verify from expression ( fllf ) that for d > 2 
and n > 1 we have da2n+i(p) I dp > 0. Negative val- 
ues of <72n+i(0) may thus become positive due to p. In 
particular, for d = 3 we obtain a^(p) = (23p — 1)/10, 
which becomes positive for p > 1/23. It then follows 
that i?3 — » oo as Pe — > oo, the footprint of the persis- 
tence of the small-scale anisotropy. 

Some remarks on the limit of large space dimensions 
are worth. One immediately realizes from Eq. ( |l4| ) that, 
for d — + oo the scaling exponents reduce to C^' q — 
—^nqdp, that means Co ' 9 = n Co' 9 = 'zCo ' 1 ' an< i thus 



the vanishing of intermittency for d — > oo. Neverthe- 
less, in this limit Co' 1 — C2' 1 an d, at the dissipative scale 
T], CT2„+i(p) = pd(2n 2 — 1/2) > 0, two clear signatures of 
small-scale anisotropy persistence. The latter may thus 
occur also in the absence of intermittency. 

Let us now examine the possible mechanisms at the ori- 
gin of the small-scale anisotropy persistence in our prob- 
lem. They can be easily grasped in two dimensions; our 
previous results being valid for all d > 2, it is reason- 
able to expect that the mechanisms we are going to show 
hold also for higher d's. Let us start from the incom- 
pressible case assuming, without loss of generality, that 
B is oriented along the y-axis (i.e. B = e y ). As we 
shall see, all our considerations will hold, a fortiori, in 
the compressible case. 

It can be shown that the magnetic field can be repre- 
sented by the (scalar) magnetic flux function in the form 
B = dip x e z , where ip satisfies the passive scalar equa- 
tion forced by a large scale gradient G = B° x e z . One 
finds 



(17) 



where ^(t, x) = ip + G ■ x. Notice that, fort d = 2, 
expression (^) recovers the results of Ref. Q obtained 
directly for the scalar problem. 

An interesting feature recently recognized for the pas- 
sive scalar problem, for both synthetic advection by 
(^-correlated velocity Jl9| and Navier-Stokes advection 
po| p3| , is related to the formation of "cliffs", i.e. very 
steep scalar gradients within very short distances. The 
region where such gradients develop are separated by 
"plateaus" where the scalar depends smoothly on the po- 
sition (i.e. regions of gradient expulsion). 
The emergence of this ubiquitous pattern is explained by 
considering the action of the velocity derivative matrix 
d a vp. As emphasized in Refs. [ p0|j2l] ], scalar gradients are 
weak in the elliptic regions of the velocity field (in two 
dimensions they correspond to purely imaginary eigen- 
values of the velocity derivative matrix) where its 'rota- 
tional' character inhibits the formation of strong scalar 
gradients. On the contrary, in the hyperbolic regions the 
flow is almost 'irrotational': the alignment of scalar gra- 
dients with the direction of the eigenvector corresponding 
to the most negative eigenvalue of the velocity derivative 
matrix (roughly, the direction corresponding to compres- 
sion along one direction) is not discouraged, and actually 
observed p5|, and strong scalar gradients develop. 



When the scalar is forced by isotropic injection mech- 
anisms, no preferential direction arises and the intense 
gradients are randomly oriented. The final result is a 
small scale isotropic statistics. The situation changes 
in the case of non-isotropic injection mechanisms like 
the one encountered here. In this case strong gradi- 
ents are oriented along G and, as very recently shown 
in Refs. [^2|,^3|, small-scale isotropy is consequently not 
restored for the scalar field. Exploiting the relation be- 
tween the magnetic field and the magnetic flux function, 
we can conclude that extreme magnetic fluctuations have 
a tendency to occur preferentially along the direction B , 
the origin of the observed small-scale anisotropy. 

In the compressible case both eigenvalues can be neg- 
ative Compression may thus occur in both directions, 
enhancing the formation of fronts in the magnetic flux 
function. 

In conclusion, we presented a simple model of MHD 
turbulence where, by varying the degree of compressibil- 
ity of the velocity field, the persistence of anisotropy is 
detectable both from the hierarchy of inertial range ex- 
ponents and from the divergence of skewness factor with 
the Peclet number. Although our results were obtained 
on the base of a specific model, they seem to be rather 
general: a similar behavior is also observed for a passive 
scalar advected by a compressible velocity field H (how- 
ever, the flattening is less pronounced for the latter). 

The results presented here give the first evidence of a 
turbulent system displaying such twofold behavior. 
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